Parity Effect in Ground State Energies of Ultrasmall Superconducting Grains 
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We study the superconductivity in small grains in the regime when the quantum level spacing 
5s is comparable to the gap A. As Se is increased, the system crosses over from superconducting to 
normal state. This crossover is studied by calculating the dependence of the ground state energy of 
a grain on the parity of the number of electrons. The states with odd numbers of particles carry an 
additional energy Ap, which shows non-monotonic dependence on Se. Our predictions can be tested 
experimentally by studying the parity-induced alternation of Coulomb blockade peak spacings in 
grains of different sizes. 



PACS numbers: 73.23.Hk, 74.20.Fg, 74.80. Bj 

The standard BCS theory § gives a good description 
of the phenomenon of superconductivity in large samples. 
However, it was noticed by Anderson in 1959 that as 
the size of a superconductor becomes smaller, and the 
quantum level spacing in the sample Se approaches the 
superconducting gap A, the BCS theory fails. The in- 
terest to the superconductivity in such ultrasmall grains 
was renewed by recent experiments by Ralph, Black and 
Tinkham HQ, who fabricated and studied nanometer- 
scale aluminum grains. In qualitative agreement with 
the prediction (2) , they demonstrated B the existence of 
superconducting gap in relatively large grains, with es- 
timated level spacings Se w 0.02 and 0.08 meV smaller 
than the superconducting gap A w 0.31 meV, whereas no 
signs of superconductivity were observed || in smaller 
grains, Se « 0.7 meV. These experiments raise a theo- 
retical question about the nature of the crossover from 
superconducting to normal state in ultrasmall particles 
with level spacings Se ~ A. 

This problem was addressed in two recent theoretical 
papers. J. von Delft et al. M explored the BCS gap 
equation in a finite-size system with equidistant discrete 
energy levels, and found that as the level spacing is in- 
creased, the superconducting gap of the grain vanishes 
at a certain critical value of Se, which is of order A and 
depends on the parity of the total number of electrons 
in the grain. Smith and Ambegaokar || extended the 
treatment of Ref. || to take into account Wigner-Dyson 
fluctuations of the energy levels in the grain. 

It is worth noting that the theories ||,|| treat the su- 
perconductivity in small grains within the selfconsistent 
mean field approximation for the superconducting order 
parameter. Although this approximation works well for 
large systems, one should expect the quantum fluctua- 
tions of the order parameter to grow when the level spac- 
ing Se reaches A. In this paper we present a theory of 
superconductivity in ultrasmall grains which includes the 
effects of quantum fluctuations of the order parameter. 



We show that the corrections to the mean field results 
which are small in large grains, Se <SC A, become impor- 
tant in the opposite limit, Se ^S> A. 

The superconducting gap A studied in Refs. d](| is 
not well defined in the presence of quantum fluctuations. 
Therefore, we must first identify an observable physical 
quantity which characterizes the superconducting prop- 
erties of small grains. The most convenient such quantity 
for our purposes is the ground state energy of the grain 
Em as a function of the number of electrons N. More pre- 
cisely, we study the so-called parity effect in ultrasmall 
grains, which is described quantitatively by parameter 



Ap — E 2 i+i — - [E21 
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In the ground state of a large superconducting grain with 
an odd number of electrons, one electron is unpaired 
and carries an additional energy Ap = A. This re- 
sult is well known in nuclear physics and was recently 
discussed in connection to superconducting grains in 
Refs. The parity effect was demonstrated exper- 

imentally in Refs. ^,0, where the Coulomb blockade 
phenomenon ifllf in a superconducting grain was stud- 
ied. In a Coulomb blockade experiment one can measure 
Ap explicitly as the difference of spacings between three 
neighboring peaks of linear conductance through a su- 
perconducting grain. 

We describe the grain by the following Hamiltonian: 



H 
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Here k is an integer numbering the single particle energy 
levels £fc, the average level spacing (e^+i — £fe) = Se, op- 
erator etfc CT annihilates an electron in state k with spin a, 
and g is the interaction constant. In Eq. (^|) we assume 
zero magnetic field, so that the electron states can be 
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chosen to be invariant under the time reversal transfor- 
mation EJ], We include in Eq. (0) only the matrix ele- 
ments of the interaction Hamiltonian responsible for the 
superconductivity; the contributions of the other terms 
are negligible in the weak coupling regime g/5s -C 1 we 
consider. Finally, we did not include in Eq. (||) the charg- 
ing energy responsible for the Coulomb blockade, as its 
contribution to the ground-state energy is trivial. 

In the absence of interactions, g — 0, the parity pa- 
rameter Ap can be easily calculated. Indeed, the ground 
state energy En is found by summing up N lowest single- 
particle energy levels. This results in E 2 i+\ = E 2 i + £ i±i 
and E21+2 = E 2 i + 2e; + i. Substituting this into Eq. (||), 
we find that without the interactions Ap = 0. 

For weak interactions one can start with the first-order 
perturbation theory in g. In this approximation an elec- 
tron in state k interacts only with an electron with the 
opposite spin in the same orbital state k. Thus when 
the "odd" (21 + l)-st electron is added to the grain, it 
is the only electron in the state I + 1 and does not con- 
tribute to the interaction energy, 5E 2 i+i — 5E2i- The 
next, (21 + 2)-nd electron goes to the same orbital state 
and interacts with it: 8E21+2 = SE21+1 — g. From Eq. (|l|) 
we now find 
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One should note that the result (||) is not quite satis- 
factory even in the weak coupling case g/5e <C 1. Indeed, 
the low-energy properties of a superconductor are usually 
completely described by the gap A. The interaction con- 
stant g is related to the gap A in a way which depends on 
a particular microscopic model, so the result (|^) cannot 
be directly compared with experiments. 

This problem can be resolved by considering correc- 
tions of higher orders in g, which are known ]l2] | to 
give rise to logarithmic renormalizations of g. In the 
leading-logarithm approximation the renormalized inter- 
action constant is found 11121 as 



9 = 



1 



f ' 

d£ D 
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Here Do is the high-energy cutoff of our model, which has 
the physical meaning of Debye frequency, and D <C -Do 
is the low-energy cutoff. At zero temperature, D ~ 5s. 
Taking into account the relation between the gap in 
a large grain A and microscopic interaction constant, 
A ~ Doe~ Se / 9 , we find with logarithmic accuracy g = 
Se/ \n(ds/ A). Finally, substituting the renormalized in- 
teraction constant into Eq. (||), we get 

Ap = ^ f , A«Se. (5) 

Unlike the first-order result (^), Ap is now expressed in 
terms of experimentally observable parameters A and 5s 
rather than model-dependent interaction constant g. 



It is instructive to compare Eq. (||) with the results of 
Refs. |,|§. In a very small grain with 5s 3> A, the mean 
field gap studied in Refs. [|[D vanishes, and no parity ef- 
fect is expected. On the contrary, our result (^|) predicts 
that in small grains the parity effect is stronger than in 
the large ones. This behavior is due to the strong quan- 
tum fluctuations of the order parameter which persist 
even when its mean field value studied in Refs. |||]| van- 
ishes. The physics of the fluctuations of the order param- 
eter is hidden in the renormalization procedure leading 
to Eq. (^). Below we present a different technique, which 
explicitly shows the role of the fluctuations. It will allow 
us to rigorously derive Eq. (@) and to study the fluctua- 
tion corrections in the case of large grains, 5s <C A. 

A convenient way to treat the fluctuations of the or- 
der parameter is by using path integral technique [ffsj . 
This approach gives an exact expression for the grand 
partition function of a superconductor: 



^,T)=Trexp(-^^V # = X>L< 
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Here fi is the chemical potential, T is the temperature, 
and N is the operator of the number of electrons in the 
grain. At T — > the dominating term in Z(fi, T) cor- 
responds to the ground state of the grain with a certain 
number of electrons: 

Z([i, T -f 0) = e - n( " )/T , = mm{E N - fiN}. (7) 

Thus we can find the ground state energy Em by studying 
the grand partition function (||). 

One problem with this method of calculating En is 
that because of the parity effect ([!]) with Ap > the 
odd charge states do not contribute to Z(/j,,T — * 0). To 
find E21+1 let us consider the effect of interactions on the 
unperturbed ground state of 21 + 1 electrons. Since the 
state l + l is filled with one electron, the interaction term 
in the Hamiltonian (Q) can neither create nor destroy a 
pair in this state. Thus E21+1 can be found as 



E21+1 — £1+1 + E 2 



(8) 



where E21 is the ground state energy of a grain with 21 
electrons for the system (^) with state k = l + l excluded. 

The idea of the path integral approach Jll| is to re- 
place the formulation (j|) of the problem in terms of 
electronic operators a^ a by an equivalent formulation in 
terms of the superconducting order parameter A(r). The 
latter is introduced as an auxiliary field for a Hubbard- 
Stratonovich transformation splitting the quartic inter- 
action term in Eq. (^) into quadratic pair creation and 
annihilation operators. Then the trace over the fcrmionic 
variables can be calculated, and one finds 



(9) 
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where the action S[A] is defined as 

1 fVT 



^[A] = -E 



TrlnGr 1 - — 
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A(r)\ 2 dT. (10) 



Here ^ — e k — ^, and the inverse Green's function 



x6(t-t'), 



A(t)(7 h 



A*( 



(11) 



where a ± = a x ± ia v , and a x ' y ' z are the standard Pauli 
matrices. G -1 satisfies antiperiodic boundary conditions: 
Gfc 1 (r + T- 1 ) = -G,T 1 (r). 

Unlike in the case of large superconductors , the or- 
der parameter A in Eqs. (p|)-(pl]) does not depend on the 
coordinates, and thus the contributions of different states 
k in the action (|l0|) decouple. This results from the sim- 
plified form of the interaction term in the Hamiltonian 
(^) . The space fluctuations of A are negligible for grains 
smaller than the coherence length of the superconductor; 
this condition is well satisfied in ultrasmall grains. On 
the other hand, the time fluctuations of A accounted for 
in Eqs. (0)-(O) lead to the corrections to the mean field 
BCS theory and are studied below. 

First we consider the regime of weak interactions, 
A <C Se. In this case the A-dependent terms can be 
considered to be a small perturbation V — Acr + + A*a~ , 
and one can formally expand the action (|l^) in power 
series in V using 



oo 

Trln(G 1 - V) = TrlnG^ 1 - ^ -Ti(G Q V) 

i=i 3 



(12) 



The first-order term vanishes because matrix V is off- 
diagonal, so we study the quadratic in A contribution 
to the action. The calculations are more convenient to 
perform in terms of the Fourier components A m of the 
order parameter, defined in a usual way: 



A(r) =T^A m e 



— tUJ m T 



— 27rTm. 



(13) 



The calculation of the second-order contribution to the 
action (n(t) is straightforward and gives 



SS = Tj2 1 -^^\A m \^ 
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a(E) = g 



V 2 ^ 



E 



(14) 



The functional integral (^) is now easily evaluated by 
integrating over the real and imaginary parts of each A TO . 
We normalize the result for the partition function Z by its 
value Zq for non-interacting system, which corresponds 
to a = in Eq. Q, 



Z(H,T) 



n- 

m 

nn 



i 



■ a(iuj m ) 
2|fc - iw„ 



n 



sinhfo/T) 
sinh(4/T)' 



(15) 



Here £fc are defined by 1 — a(2£fc) = 0. Assuming weak 
interactions, A <C Se, we find £k = £,k + S£k, where 



sgn^fe 



2 1 - 2V„ 



(16) 



We can now compare the result (|l~5|) with the definition 
(0) of tl(fi) and find 



(17) 



One can easily see that for e; < /j, < this correction 
does not depend on ^i. According to (0) such (5fi should 
be interpreted as the correction 8E21 to the ground state 
energy of 21 electrons present in the grain in this range 
of fi. We then use the rule (^) to find SE21+1. To do this 
we exclude the state k = 1+1 from the sums in Eq. jl^ ) 
and ([l?]) and calculate SE' 2l . The result for Ap coincides 
with Eq. (|). 

We now turn to the case of stronger interactions, 
i > fe. A good starting point in this regime is the 
standard BCS theory which corresponds to a mean- 
field approximation for the order parameter in the path 
integral approach [13|. Substituting a time- independent 
A into the action (10), one finds 



n(^) = J2(^-e k ) + -\A\ 



(18) 



Here e k = ((% + | Ap) 1 / 2 , and the value of |A| must be 
chosen in a way which minimizes fl. This means that | A| 
is determined from the usual BCS equation: 



E — = -• 



(19) 



In the continuous limit 5e/A — > one can apply the rule 
(||) and find that the exclusion of one state from the sum 
over k in Eq. ( |T§| ) results in the energies of odd charge 
states exceeding those of the even ones by Ap = A. 

To find the corrections to Ap due to the finite Se a 
more careful treatment of the mean-field approximation 
is required. One can easily see that not only a time- 
independent Ao(r) = |A|, but also any path Am(t) — 
\A\e l2nMTT with integer M is a minimum of the action 
( |To| ) which must be taken into account. It is convenient to 
treat the path Am (t) in Eq. ( pd| ) by performing a gauge 
transformation U = exp(iirMTTa z ), which eliminates 
the time dependence of Am(t) and shifts the chemical 
potential fx -> fx-iirMT @. Thus instead of Z = e~°/ T 
the partition function at T — > is now 
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Z(jjl,T) =J2e~ n{ "~™ MTyT 

M 



M 



where fl(fi) is given by Eq. (|18|). We have expanded 
n(/j, — iirMT) in Taylor series in iirMT, and neglected 
the terms vanishing at T — > 0. It is now obvious that 
Z(/j,, T) = unless the derivative of is an even integer: 



-21. 



(20) 



Thus the mean field approximation can be applied only 
for discrete values of chemical potential /i 2 / correspond- 
ing to solutions with 21 electrons. For odd number of 
electrons H21+1 is found as H21 in a system with one state 
k = I + 1 at the Fermi level excluded. At A ^> Se one 
always gets ^n+i — UN — Se/2. 

To find Ap we substitute in Eq. (Q) the ground state 
energy as En = f2(/ijv) + I^nN . The contribution of the 
second term to Ap is — Se/2. In evaluating the contri- 
bution of f2(/Azv) one has to take into account the de- 
pendence of [i jy on N and the suppression |?],|| of the 
self-consistent gap in Eq. ( |l8| ) A Q dd = A — Se/2 for odd 
N due to the exclusion of one state k from the gap equa- 
tion (|l9|). Combining all the contributions, we get the 
following mean- field result: 



Se 



Se < A. 



(21) 



It is interesting that a similar quantity Ap = — E21 + 
(i?2i-i + E2i+i)/2 is unaffected by finite level spacing up 
to the terms linear in Se, i.e., Ap = A. Thus at Se <C A 
we have Ap = A Q dd and Ap = A even . 

Though the result Ap = A can be obtained from the 
mean field theory, an evaluation of corrections to it due 
to the level spacing requires taking into account the ef- 
fects of the fluctuations of the order parameter. One can 
find the contribution of the fluctuations by expanding the 
action near A(r) = |A| using Eq. JI2]). The second-order 
correction to the action is [[LSI 



SS = Tj2\s^U 2 + i6 I m^ 1/2 \ E 



2e fe (w m 7 m - 2i£ fc ) ' 



where S^ 1 are the Fourier components of the real and 
imaginary parts of the fluctuation A(t) — |A|, and j m = 
(1 + 4|A| 2 /o;^) 1 / 2 . Now we evaluate the path integral 
(||) by integrating over all S^ 1 and find the contribution 
of the fluctuations to £!(/i), 



Sn(fi) 



du> 
2?r 



In 



E 
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2e fc (Vw 5 + 4lA]*-2i& 



A comparison of the expressions for <5f2 in the cases of 
even and odd numbers of electrons shows that they co- 
incide up to the terms of order Se. Thus the fluctuations 
of the order parameter do not affect the result (|2l]) . 

Finally, we discuss the mesoscopic fluctuations of Ap 
given by our results (j|) and (^T|). It is clear from 



Eq. ( |l6| ) that unlike Se in the numerator of Eq. (||), 
the one in the argument of the logarithm is sensitive 
to the Wigner-Dyson fluctuations of Thus the rel- 
ative mesoscopic fluctuation of the result (||) is small, 
SAp/Ap ~ l/ln(fc/A). It is also interesting to com- 
pare the mesoscopic fluctuation of the gap A originating 
from the level fluctuations in Eq. ( [l9| ) with the small cor- 
rections in Eq. (21). One can easily express [lq] the cor- 



rection to A in terms of the correction to the density of 
states z^(£) in Eq. (|l9|). Then the mean-square fluctuation 
of the gap is found using the well-known results for the 
correlator {v(g)v (€')), and we get ^((<5A) 2 ) = Se/ny/2. 

In conclusion, we have studied the parity effect (jl]) in 
the ground state energies of an ultrasmall superconduct- 
ing grain. Although the quantum fluctuations of the or- 
der parameter can be neglected for large grains, Eq. (|2l]), 
they play a crucial role in small grains, Eq. (||). As the 
size of the grain decreases, the parity effect first weak- 
ens, Eq. (^l|), but then starts increasing, Eq. (||). Thus 
we expect a minimum of Ap at a certain size of the grain 
such that Se ~ A. 
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